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We analyze the possibihty of the nanoscale phase separation manifesting itself in the formation of 
ferromagnetic (FM) polarons (FM droplets) in the general situation of doped anisotropic three- and 
two-dimensional antiferromagnets. In these cases, we calculate the shape of the most energetically 
favorable droplets. We show that the binding energy and the volume of a FM droplet in the 
three-dimensional (3D) case depend only upon two universal parameters J — (J^ + Jy + Jz)S^ 



and teff = {txtyt^Y^^ , where J and t^ff are effective antiferromagnetic (AFM) exchange and 
hopping integrals, respectively. In the two-dimensional (2D) case, these parameters have the form 
J = {Jx + Jy)S^ and t^ff = {tj:tyy^^. The most favorable shape of a ferromagnetic droplet 
corresponds to an ellipse in the 2D case and to an ellipsoid in the 3D case. 
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I. INTRODUCTION 

The problem of electronic phase separation with the 
formation of ferromagnetic (FM) or paramagnetic (PM) 
spin polarons (magnetic droplets or ferrons) due to the 
self-trapping of charge carriers in an antiferromagnetic 
(AFM) matrix became very popular nowadays, espe- 
cially in the studies of high-Tc superconductors and the 
systems with the colossal magnetoresistance (such as 
LaMnOa manganites doped by Ca). For isotropic ma- 
terials, the size and shape of FM droplets was evalu- 
ated in several papers beginning from the seminal work 
of Nagae\ii, for more detail see Ref. ^ The charac- 
teristic size of a FM droplet turns out to be of the 
order of 15-20Aand its optimum shape in isotropic 
3D manganites is a spherical one. Later on, Kagan 
and Kugel'^ analyzed the case of layered manganites 
(like(La,Ca)„+iMn„03„+i) and demonstrated that the 
droplets with the lowest energy have the ellipsoidal 
shape. The FM droplets of cylindrical shape considered 
first by Nagae"\i^ for this class of manganites correspond 
to a higher energy. Currently, the phase separation in 
anisotropic materials was also addressed in connection 
with low-dimensional organic compounds^ and quasi- 
one-dimcnsional magnets such as BaCoOj^^. Magnetic 
polarons in doped one-dimensional (ID) AFM magnetic 
chains were recently considered in Refs. I^lgl liol where 
the possibility of rather long-range magnetic distortions 
around the polaron was demonstrated. 

Another possibility to have a strongly anisotropic sit- 
uation arises when we take into account an interplay be- 
tween the microscopic phase separation and charge or- 
dering (stripe formation), include the Jahn- Teller type 



of effects (orbital degrees of freedom), or consider stable 
crystallographic distortions. In these cases, the quasi- 
ID zig-zag or ladder structures are often observed in the 
corresponding systemsii*i^. 

In this paper, we present the calculations concerning 
the shape and size of FM droplets in anisotropic two- 
dimensional (2D) or three-dimensional (3D) cases when, 
generally speaking, the electron hopping integrals t^, ty, 
and tz along x, y, and z directions, as well as the AFM 
exchange integrals Jx, Jy, and J^ are different. We get 
that, by analogy with the situation in layered mangan- 
ites^, the most favorable shape of a FM droplet is an 
ellipsoidal one. Moreover, the binding energy and the ef- 
fective volume of the droplet are expressed only in terms 
of universal averaged parameters J = {Jx + Jy + Jz) S'^ 

1/3 

and tf.ff — {txtytz) ■ These results are interesting, in 
particular, in relation to the neutron scattering experi- 
ments giving an indication of the existence of FM clusters 
with different shapes in perovskite and layered mangan- 

Our paper is organized as follows. First, we consider 
the purely 2D situation and find the most favorable shape 
of a 2D ferron comparing the energies of elliptical and 
rectangular droplets in the general anisotropic 2D case: 
tx ^ ty and Jx ^ Jy We find that in two dimensions, 
the minimal energy corresponds to the elliptical shape. 
Then, we include the third dimension (J^ and tz) and 
compare the energies of the cylinder and ellipse in the 
case when both of them have the optimum elliptical shape 
of the 2D cross-sections. We find again that the minimal 
energy in the 3D case corresponds to the ellipsoidal shape 
of FM droplets. At the end of the paper, we provide some 
discussions and conclusions. 
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Figure 1: 2D anisotropic doped antiferromagnet with the 
square lattice. Ja; 7^ Jy are AFM exchange integrals, t^^ 7^ 
are the electron hopping integrals, and Ly define the 
volume of the rectangular ferron in 2D. 



II. THE SHAPE OF FM DROPLETS IN THE 
ANISOTROPIC 2D CASE 



A. A rectangular ferron 

Let us first consider a rectangular FM droplet (ferron) 
located at the square lattice with the intersite distance 
a. Its characteristic sizes along x and y axes are and 
Lj,, respectively. The dimensionless volume of such 
a ferron can be defined as SI = LxLy/a?. The kinetic 
energy of charge carriers (electrons or holes) within the 
FM droplet is 



Ekin ^ -2txn - 2tyn + eon, 



(2) 



where n is the concentration of charge carriers and Sq is 
a binding energy corresponding to the first (the lowest) 
level in the rectangular potential well. The latter can be 
found by solving the corresponding Schrodinger equation 
(see Refi) 



where 



tlkin — a \tx ^_^2 + h Qy2 



(3) 
(4) 



For the case of a well-defined ferron (without an ex- 
tended tail of magnetic distortions, which we will con- 
sider throughout the present paper, the corresponding 
boundary conditions have the form 



Let us first consider the anisotropic 2D case. In this 
case, there are two different electron hopping integrals 
tx 7^ ty and two different constants of the AFM exchange 
interaction ^ Jy, see Fig. ^ To some extent, this 
case has a lot of similarities with two-leg ladder systems 
rather popular nowadays (see, for example, Ref. [I^ and 
references therein). 

Throughout this paper, we consider the Kondo-lattice 
model with the Hamiltonian 

H ^ Jh^^ Sia^ + ^ Ja^iSj + ^ taclcj, (1) 

« (U>a {ij)a 

where, c| and Ci are electron creation and annihilation 
operators at site i, a = {x,y} for a square lattice in 2D, 
{ij)a denote the neighboring sites in the lattice along 
the a direction, cr^ = c|crci is the spin of a conduction 
electron (ct is the Pauli matrix). Si is a local spin, Jq, are 
AFM exchange integrals, ta are the hopping integrals for 
conduction electrons, and parameter Jh corresponds to 
the Hund's rule coupling between a local spin S and a 
spin of a conduction electron. 

We work in the double-exchange limit, which implies 
that Jh ^ {tx,ty} ^ {Jx, Jy}- In this case, the ground 
state of the system is unstable toward the nanoscale 
phase separatioijiiSi^ with the formation of FM polarons 
inside the AFM matrix. Let us now evaluate the total 
energy of the phase-separated state for different shapes 
of ferrous possible in the 2D case. 



■i:{x^L^,y)^^{x,y^Ly)^Q. (5) 

Hence, 

, . . TTX . Try 
^>(x,y) = sm — sm— (6) 

Ly 

and 

/7ra\ ^ /7ra\ ^ , , 

''-'At:) ^'•{tJ ■ 

Now, we can pass to the evaluation of the potential 
energy given by the terms related to the AFM exchange 
interaction. In the domains with the ferromagnetic or- 
der (ferrous), the AFM exchange leads to the positive 
contribution to the total energy 

Ep,ti = 2{J^ + Jy)S\^^. (8) 

For the AFM regions, which are free of ferrons, the cor- 
responding contribution to the energy can be written as 

Ej,,,,^~2{J, + Jy)S''(^l-n^y (9) 
Hence, the total potential energy yields 
Epot = -2 (J, + Jy) ^2 + 4 (J, + Jy) S^n^^. (10) 
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As a result, the total energy related to the formation 
of FM droplets has the following form 

Etot = -2 [t^n + tyn + ( J^, + Jy) 5^] + 



2 / \ 2 



Lx 



The minimization of energy with respect to 
and Lj, gives 



dEt. 



dEtot 



~2nt 



+ 4(J, + Jj,)n52% ^0. (12) 



A solution to equations (|12(l reads 



L 

txT^ =2 ( + Jy) 5*^ 



l4 



( Jr + Jy) S 



(13) 



LyLx 



■n trt. 



x''y 



Multiplying two equations (fT!^ by each other, we find 

(14) 

(15) 



a- J A{Jx + JyYS^ 

Now, introducing notation 



te// — {txty) , J — ( Jx + Jy) S 



we find 



LyLx 



4J2 



Thus, the dimensionless volume (area) f7 of a 2D ferron 
can be written as 



V2 \ J 



1/2 



(16) 



We get quite a remarkable relationship expressing the 
volume of a 2D ferron in terms of the teff/J ratio. 

Correspondingly, the minimized total energy lll|l takes 
the form 

Etot = -2 (txTi + tyji + j) + 47rV2n (teff.T)^^'^ ■ (17) 

Introducing the energy of FM polaron by the relation- 
ship 



Epoi = Etot + 2 {txTi + tyn + J) 
we get finally 



Epoi = SnilJ = 47rV2n {teffJ 



1/2 



(18) 



(19) 



It is again worth to notice that the energy of the FM 
polaron in 2D depends only upon the product of t^ff 
and J. 



B. An elliptical ferron 

Now, we can consider the energy a two-dimensional 
FM polaron having the shape of an ellipse. For the same 
characteristic sizes (principal axes)of the ferron, its vol- 
ume in the case of an ellipse is f2 = TiLxLy/ a? . The 
corresponding kinetic energy is again given by Eqs. 
(0J. To solve the Schrodinger equation in this geometry, 
we should transform an ellipse to a circle. This could 
be done, for example, by the dilatation along the y axis: 
y = ijy/ ty/tx- Then, we have 



Hhin 



where and R 



-a^tx 



9y2 



-a t,rAs 



(20) 



r2 + 2;2andA^ = ^ + i^isthe 
radial part of the Laplacian operator in 2D. Thus, the 
ellipse /Li. + y"^ /L"^ = 1 in the 'old' x,y coordinates 

transforms to the + y^ — Rmax circle in the 'new' 
X, y coordinates. From equation for the circle in terms of 
'new' coordinates x, y, it is clear that Rmax = Lx- Hence, 
we have 



Lx — Ly \l — — Rr, 



(21) 



and the ferron volume in the initial ('old') coordinates 
reads 



LxLy 

il = TT — 2- 



0? V 



(22) 



In this solution to the Schrodinger equation (PJ 

has the form \1/ = Jo(fci?), where Jo is the Bessel function 
of zeroth order. The boundary condition Jo{kRmax) — 
yields kRmax — jo,i = 2.404 « 37r/4, where jo,i is the 
first zero of function Jq . This means that 



eo = txtt^k^ — t. 



Then, we have 



Jo.iQ 

Rmax 



= tx 



Lx 



(23) 



Etot = -2 [txTi + tyu + {Jx + Jy) S^] + Epoi, (24) 



where 



Epoi^n[^-^ytx + ^ [Jx + Jy) S'n^ ^[f. (25) 



Lx 



a V 



The minimization of polaron energy ()25|l with respect to 
Lx gives 



dEpoi 
dLx 



n , i0,l"^ , 7 I 7 \ C2 jty 

-2ntx — — \-8{Jx + Jy)o n- ' 



LI 



V tx 



= 0. 

(26) 
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Thus, we have (see Eqs. and (PSjl l 



f 



""Jo, 1^2; 



(27) 



Introducing again t^ff and J defined by Eq. (|15|l . we 
find 



jo.lV^ f tefj_ 

J 



1/2 



(28) 



So, the ferron volume is again expressed as a function of 
the universal ratio te///J. Comparing expressions 
and H28|l for the volumes of rectangular and elliptical fer- 
rons, we find 



ellipse 



Jo,i 



rectangle 



0.96 < 1. 



(29) 



This means that the elliptical ferron is a more compact 
object (i.e. it has a smaller volume) in comparison to the 
rectangular ferron. Accordingly, the energy of elliptical 
magnetic polaron can be written in the following form 



1/2 



(30) 



We can see again that the ferron energy depends only 
upon the product of te// and J. Finally, we can com- 
pare the ferron energies for the cases of rectangular and 
elliptical shapes using Eqs. ((T^ and l|Sn|) 



e: 



^ellipse 
pol 



Jo,i 



ellipse 



E 



'rectangle 
pol 



2tt ^ 



0.96 < 1 



rectangle 



(31) 



We see that the ratio of energies turns out to be iden- 
tical to the ratio of the volumes. Thus, the elliptical 
shape of the ferron is more favorable in energy than the 
rectangular shape. First of all, this is caused by a more 
compact structure of the elliptical ferron. Another cause 
was emphasized in Ref. (sl The thing is that the ellip- 
tical shape of ferron in 2D has a close similarity to the 
one-electron spectrum characteristic of the empty square 
lattice: Sp — p^/2mx +Py/2'my, where mx/2 = txO^ and 
my/2 — tyO^. Concluding this section, we can say that 
the elliptical shape is the shape most favorable in energy 
for a FM droplet in doped anisotropic antiferromagnets 
with the 2D square lattice. 



III. THE SHAPE OF FM DROPLETS IN THE 
ANISOTROPIC 3D CASE 

Now, we can include the third dimension (which means 
the inclusion of Jz and tz) and consider the shape of a 
FM droplet in a doped anisotropic antiferromagnet with 
the 3D cubic lattice. Of course (having in mind the re- 
sults of the previous section), we have to consider FM 
droplets with the 2D cross-section most favorable in en- 
ergy. In other words, we consider the 3D droplets hav- 
ing the shape of an ellipse in the x, y plane. Then the 



problem effectively reduces to the comparison of energies 
and volumes of a cylinder and of an ellipsoid of rotation 
both having an elliptical cross-section with dimensions 

= LyJ^ (see Eq. (jSB). 



A. FM droplets of cylindrical shape 

First, let us consider 3D FM droplets of cylindrical 
shape. The volume of such a droplet can be written as 



^ = 1 — - 2 \i , 

In this case, the total energy has the form 




(32) 



Etot = -2 [txTl + tyU + tzU + [Jx + Jy + Jz) S"^] + 



The polaron energy 

Epol = Etot + 2 [txTl + tytl + tzTL + {Jx + Jy + Jz) S"^] 

(34) 

is given by the expression 



E, 



pol — ''X 



4(Jx + Jy + Jz)S^mT-^^ /^^-^^ 



(35) 



The minimization of polaron energy H35() with respect 
to Lx and yields 



9 Epol „, Jo,i^ 
— ^— = —2txn — — — 

^x 

-2tzn- 



dLx 



dEpoi 



dLz 



LI 



.8Jn^%^^./^ = 0, 
a a \ tx 

AJnir^.f^^O. (36) 



where we introduced the effective exchange integral for 
the 3D case 



J^iJx + Jy + Jz)S^, 

From Eqs. we get 



(37) 




(38) 



Squaring the second equation in H38|) and dividing the 
result by the first equation, we exclude and obtain 
the following expression for Lz 




1/5 



(39) 
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Substituting Eq. (|39|l to the first equation in l|38l) , we get 



Hence the energy of the elliptical FM polaron can be 
written as 



1/5 



(40) 



Hence, we have 



1/5 



V2 1 ttV 



(41) 



Using Eqs. (|39|l and (|41|l . we find volume ^cyi of the 
cylindrical fcrron 



(42) 



\l tx 2 V 

Introducing the effective hopping integral for the 3D case 



we can rewrite Eq. H42|l as 



cyl 



(W)'^' ftef_f\ 



\ 3/5 



(43) 



(44) 



Similar to the the 2D case, we see that the ferron vol- 
ume in 3D is also a function of tejf/J ratio, where the 
effective parameters are given by Eqs. H37(l and (|43|l . 

Substituting expressions (jSHJ, , and for L^, 

Lx , and U,cyi , respectively, to the energy of a FM polaron 
we get 



Ep,i = lOnm = 5n(^jo,i)^/' i^lffJ^f^^ 



(45) 



We see that the polaron energy in the 3D case again 
depends on the universal parameters t^ff and J, but 
the specific form of this dependence is slightly different: 

/ -\l/5 -1/2 

(tg^jJ^j in 3D as compared to ifeffJ) in 2D. 



B. FM droplets of ellipsoidal shape 

Here, we calculate the volume and the energy of the 
FM droplet having the ellipsoidal shape. The volume of 
the ellipsoidal droplet in the 3D case is 



4 LxLyLz 



(46) 



In this case, the total energy of the system has the 
form 



Etot = -2 [txTl + tytl + tyU + [Jx + Jy + Jz) S ] 

+eon + AiJx + Jy + Jz) S'^iln. (47) 



Epoi — Son + AJfln, 



(48) 



where we again introduce J defined by Eq. H37(l . 

As above, energy £o can be found by solving the cor- 
responding Schrodinger equation 



(49) 



where 



/ ()2 q2 q2 \ 

Using the the dilatation along the y and z axes: y = 
V\fuji'x and z = z^tzjtx, we get 

ilkin = -txO^Aj^ (51) 

in the'new' coordinates x, y, and z. Here, we have — 



' y^ + 2? and A5 = 7^ + 2i^ is the radial part of 



d 



1 a 



R dB. 



the Laplacian operator in 3D. 

In these coordinates, a droplet is confined within a 
sphere of radius Rmax = Ex- Accordingly, we have 



Ly\ ix /iy — Lz \/ tx/tz — Lx — R-n 



(52) 



and the droplet volume expressed in terms of initial 
('old') coordinates reads 



4 LxLyLz 
^ieii — t:'"' 5 — 



4 

— vr 

3 \ a 



Ex \ (tytz ) 



1/2 



tx 



(53) 



A solution to the Schrodinger equation (|49|l has the 
form 



*(fci?) 



sm{kR) 
(kR) 



The boundary condition ^>{kRr, 
kRmax — Ti"- Hence, we find 

^2 

So ~ txO^k^ = txO^ — 



(54) 

„) = yields 
txa^-j (55) 



R2 -^"^ ^2 

^max ^ 

and the energy of ellipsoidal FM polaron takes the form 
Epol = tx—p^n + Un-TT [ — ) ^^-f— (56) 



ix 



The minimization of polaron energy (|48|l with respect to 
Lx yields 



dLx 



= -2t, 



^^n+mj.n(^\ iV^Zi^O (57) 
Li \ a a tx 



As a result, we get the expression for Lx 



T , TT t 

Ex ^ a \ — 



J y^tytz ) 



1/5 



(58) 
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Substituting Eq. lf55|) to itS^ . we find the volume of an 
ellipsoidal droplet 



^ell = 



J 



(59) 



where t^ff is again given by Eq. H43|l . We see that the 
volume of the ellipsoidal droplet is also determined by the 
dimensionless universal ratio teff/J- Dividing Eq. l(5^ 
by Eq. H44|) , we obtain the ratio of volumes for ellipsoidal 
and cylindrical droplets 



n 



ell 



^cyl 3 \2jQ 



4/5 



0.95 < 1. 



(60) 



We see that the ellipsoidal FM droplet is a more compact 
object (with a smaller volume) than the cylindrical one. 

Substituting expression (|58|l for to the polaron en- 
ergy (|56|l . we find 

8/591/5 _ , 
Epoi - lOnjn^u = 1071 {tlffP) ^ ■ (61) 

Hence the ratio of energies corresponding to two different 
shapes of ferrous is again identically equal to the ratio of 
their volumes 



^pol _ ^ell _ 4 / TT 



E2l ^cyl 3 V2jo,i 



4/5 



0.95 < 1. 



(62) 



Thus, in the 3D case, the ellipsoidal droplet has the low- 
est energy in agreement with the results of Ref. 0. 



IV. CONCLUSIONS AND DISCUSSION 

We considered the formation and the shape of droplets 
in the most general cases of doped anisotropic 2D and 3D 



antiferromagnets with arbitrary values of the electron 
hopping integrals and the AFM exchange integrals 
Ja- We found that in the anisotropic 2D case (when 
a — {x,y} and ^ tj,, Jx 7^ Jy), the most energeti- 
cally favorable shape of FM droplets is an ellipse. In the 
anisotropic 3D case (when a — {x,y,z} and we include 
into consideration the third dimension with tz and Jz), 
the most energetically favorable shape of FM droplets is 
an ellipsoid. Moreover, both the binding energy and the 
volume of FM droplets depend in the 2D as well as in the 

3D cases upon only two universal parameters t^f / and J. 

1/2 

In the 2D case, these parameters are te// = {txty) 
and J = {Jx + Jy) , whereas, in the 3D case, the cor- 
responding expressions have the form te// = {txtytz)^^^ 
and J = {J^ + Jy + Jz) 5"^. 

Note that in the present paper, we considered only the 
case of 'free' ferrous, which are not strongly localized at 
donor impurities. The study of strongly localized fer- 
rous bound to impurities, especially their shape and the 
form of the cloud of magnetic distortions related to them 
(similar to those described in Refs.lgl and flTj) will be the 
subject of a separate publication. 

Note also that the situation would be more complicated 
for FM droplets in the frustrated triangular or kagome 
lattices. This is a case, for example, in an interesting 
quasi-lD magnetic material BaCoOa, where the chains 
of Co^^ ions form a triangular latticaSii. 
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